An improved limit on the neutrino mass with CMB and redshift-dependent halo 
bias-mass relations from SDSS, DEEP2, and Lyman-Break Galaxies 
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We use measurements of luminosity-dependent galaxy bias at several different redshifts, SDSS at 
z = 0.05, DEEP2 at 2 = 1 and LBGs at z = 3.8, combined with WMAP five-year cosmic microwave 
background anisotropy data and SDSS Red Luminous Galaxy survey three-dimensional clustering 
power spectrum to put constraints on cosmological parameters. Fitting this combined dataset, we 
show that the luminosity-dependent bias data that probe the relation between halo bias and halo 
mass and its redshift evolution are very sensitive to sum of the neutrino masses: in particular we 
obtain the upper limit of < 0.28eV at the 95% confidence level for a KC DM + mv model, with 

a erg equal to erg = 0.759 ± 0.025 (Icr). When we allow the dark energy equation of state parameter 
w to vary we find w = — f .30±0.f9 for a general wCDM + mi, model with the 95% confidence level 
upper limit on the neutrino masses at "Y^rUv < 0.59eV. The constraint on the dark energy equation 
of state further improves to ui = — i.f25 ± 0.092 when using also ACBAR and supernovae Union 
data, in addition to above, with a prior on the Hubble constant from the ffubble Space Telescope. 
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I. INTRODUCTION 

Galaxy clustering at large physical or angular scales 
corresponding to the linear regime is a well-known 
probe of the cosmological parameters. The com- 
bination of cosmic microwave background (CMB) 
anisotropy data and the galaxy clustering spectrum, 
when combined with additional probes such as dis- 
tance measurements with Type la supernovae or 
baryon acoustic oscillations, are known to break var- 
ious degeneracies between cosmological parameters 
that exist when usingeither CMB data or galaxy clus- 
tering data alone [Ell]. 

While cosmological parameters are generally de- 
rived from the shape of the galaxy power spectrum by 
marginalizing over the overall uncertainty associated 
with galaxy bias that relates clustering of dark matter 
to galaxies at the linear scales, the galaxy bias, hg, it- 
self contains certain cosmological information that is 
generally ignored. The relation between galaxy bias 
and cosmology is evident in the context of the halo 
model for galaxy clustering, which provides a sim- 
ple way to relate the galaxy distribution to that of 
the dark matter halo distribution. The cosmological 
information is present in the relation between dark 
matter halo distribution and the linear density field 
and is captured by the dark matter halo bias as a 
function of the halo mass bh{M) [3|. The dark mat- 
ter halo bias contains shape information of the power 
spectrum through rms fluctuations of mass a{M). 

In the halo model, either expressed in terms of a 
halo occupation distribution (HOD; 0, H, H, 0, H, 

i, p 

(CLF; [1 



121) or a conditional luminosity function 
14,[lB|), one can relate the clustering bias 
of galaxies measured as a function of the galaxy prop- 
erty, such as the luminosity bg{L), to the bias of dark 
matter haloes as a function of the halo mass. To de- 
rive the relation between hg{L) and bh{M), it is nec- 



essary to have a proper understanding of the relation 
between an observable quantity such as the galaxy 
luminosity L and a more fundamental quantity, the 
halo mass M [H, [11, [O, Ull. This relation can be 
achieved through modeling of certain galaxy observ- 
ables, such as the luminosity function, non-linear or 1- 
halo part of the galaxy clustering spectrum, and rela- 
tions related to galaxy-mass observables from galaxy- 
galaxy weak lensing measurements [l^, [H, [2l|. In 
the context of cosmological measurements, the bg{L) 
relation as measured by the SDSS survey at low red- 
shifts has been used to constrain cosmological parame- 
ters, including the neutrino mass [H, HJ- When com- 
bined with WMAP 1-year data, SDSS power spectrum 
shape, and the SDSS hg{L) relation results in a 95% 
confidence limit on the sum of the neutrino masses of 
~ 0.54 eV. 

Beyond SDSS, several galaxy surveys that target 
galaxy populations at higher redshifts either through 
spectroscopic measurements or through Lyman drop- 
out techniques have provided measurements related to 
the galaxy luminosity functions and galaxy clustering 
power spectra or correlation functions, as a function 
of the galaxy luminosity. Among these surveys are the 
DEEP2 [2^ at z ~ 1 and Lyman-break galaxy surveys 
at z ~ 3 to 4 m, [23 . These clustering and luminosity 
function measurements can be interpreted in terms of 
a common CLF model [l^, [H, [l^ with which one 
can derive the appropriate relation to connect gal axy 
luminosity L to halo mass M at a given redshift z jl9l 
I20, 21] terms of a conditional probability distribution 
fimction P{M\L,z) [H, [l3, Hl . 

Given that clustering measurements at large linear 
scales lead to estimates of hg [L) at a redshift different 
from SDSS and we also have a mechanism to connect 
hg{L) to hh{M) through a statistical description, this 
raises the possibility of further constraining the cos- 
mological parameters than using SDSS galaxy power 
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spectrum and SDSS galaxy bias-luminosity relation 
alone. In addition to the shape information captured 
by bg{L) at each redshift, the overall evolution of bg{L) 
as a function of redshift is further sensitive to the lin- 
ear growth function of dark matter fluctuations. Since 
the growth function depends strongly on properties of 
the dark energy, such as the equation-of-state (EOS) 
relating the ratio of dark energy pressure to density, 
the combination of bg{L) measurements at several red- 
shifts raises the possibility of constraining the EOS, 
in addition to cosmological parameters that probe the 
shape of the dark matter power spectrum. 

The paper is organized as follows. In § 2 we provide 
a brief summary on how cosmological information can 
be extracted from bg{L,z) measurements by making 
use of the relations between galaxy luminosity and 
halo mass, captured by the probability distribution of 
a galaxy with a luminosity L to appear in a halo mass 
of mass M at a redshift z, P{L\M, z), from CLE mod- 
eling described in Ref. [ia,|2a|- In § 3, we describe the 
analysis of all data. In addition to bg{L) measure- 
ments at three different redshifts, we also make use 
of WMAP 5-year datajH (by updating WMAP first- 
year data used in the analysis of [1^) and shape of the 
SDSS Luminous Red Galaxy (LRG) power spectrum 
[is] at low redshifts. § 4 presents our results and we 
concluded with a summary of important constraints 
on cosmological parameters in § 5. 



II. MODELLING THE BIAS 

The clustering of bound viralized objects is biased 
with respect to that of underlying dark matter distri- 
bution and with the decreasing number density of the 
objects, the bias factor is known to increase [3]. Thus, 
bright galaxies that are in rare massive halos are ex- 
pected to be more biased that less luminous and abun- 
dant galaxies. Here, we model the relation between 
galaxy bias and halo bias following an approach sim- 
ilar to that of Ref. [2^, but using the improved halo 
bias relation from Ref. ^29^ corresponding to the el- 
lipsoidal collapse model instead of the fitting function 
for bias [i^ . We also generalize this relation to higher 
redshifts (see. Appendix A of Ref. [4]|): 



can be calculated through the relation: 



^/ab[. 



av 



bh{v{z)) = 1 + ^-[^^(a^/2)+ 
{av-^Y + b{l-c){l-c/2y' 



(1) 



in this expression 5c = 1.686 is the threshold over- 
density required for collapse of an over-density region 
and v{z) — 5c/a{M,z). The parameters a, 6 e c are 
constants and we use the values suggested in Ref. [1^ 
with a — 0.707, b — 0.5, and c = 0.6. The quantity 
a{M, z) is the rms mass fluctuation in spheres with 
radius r = {SM/iirpY^^, where M is the halo mass 
and p the mean matter density at redshift z. a{M, z) 



a'^{M,z) 



27r2 



P{k,z)W^{k)k^dk (2) 



where W{k) is the Eourier transform of ttop-hat win- 
dow function. To compute equation Q we use the 
linear matter power spectrum P(fc, z) generated by 
CAMB at redshift z for a given set of cosmological 
parameters. The dependence of bias on cosmologi- 
cal parameters is contained in the quantity v{z) = 
dc/cr{M, z) through information from the linear mat- 
ter power spectrum from a{M, z). 

As shown in several works involving measurements 
with data galaxy bias, as measured from the galaxy 
power spectrum at large physical scales corresponding 
to linear regime of clustering depends on luminosity, 
with brighter galaxies more strongly clustered than 
fainter ones [SO, HH . In the halo model, galaxies are 
expected to populate dark matter halos and these ha- 
los are already biased with respect to the density field 
bh{M,z), where M is the mass of a halo at a red- 
shift z. This is the quantity that is directly linked 
to cosmology, while bias measurements directly from 
data are a function of luminosity. We can relate the 
two thr oug h the probability distribution P{M; L, z) 
[S El) l32l | that a galaxy of luminosity L resides in a 
halo of mass M. 

If we know the P{M; L, z) then the bias at a fixed 
luminosity is given by: 



b{L,z)= / P{M;L,z)bhiM,z)dM 



(3) 



For SDSS galaxies, the conditional probability 
P{M; L, z) at low redshifts was derived based on a 
combination if SDSS galaxy luminosity function [s^ 
and luminosity-dependent galaxy correlation func- 
tions [131 that probe the non-linear, 1-halo term of 
the halo model [5| . The luminosity function is a strong 
probe of the Lc{M) relation relating the luminosity of 
central galaxies to their host dark matter halo mass, 
as well as an average scatter in that relation [Tsj , while 
the non- linear (1-halo) part of the galaxy clustering, 
either the correlation function or the power spectrum, 
establishes information related to the CLE of satellite 
galaxies. The large, linear scale clustering provides 
necessary information related to bg{L). The degen- 
eracies in the model parameters related to the CLE 
parameterization is broken with additional data such 
as the of galaxy-mass correlation function from SDSS 
galaxy-galaxy lensing measurements, similar to the 
analysis in Ref. [23j , and we make use of results from 
publicly available SD SS g alaxy-mass correlation func- 
tions [35i, [36( in Ref. In the case of DEEP2 and 
higher redshift LEG data, we again use galaxy cluster- 
ing (DEEP2: Ref. t2|; Subaru LEG: Ref. [271 and lu- 
minosity functions measurements (DEEP2: Ref. [37j : 
LEG: Ref. [26^) from the literature [1^. 
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central galaxies and the halo mass, the distributions 
have a scatter even for the central galaxy peak at the 
low-mass end. The previous analysis in Ref. [2^ ig- 
nored this scatter and described the relation between 
central galaxy luminosity and halo mass with a delta 
function and assumed simple model description with 
one free parameter to describe the same relation for 
satellites. The distributions shown in fig. [T] have addi- 
tional uncertainties due to limitations in constructing 
CLFs and when fitting to data, we allow for this uncer- 
tainty in two ways: to account for an overall system- 
atic error, we marginalize over a nuisance parameter 
b^, that scales the bias values by an overall factor and 
we include an additional error in bias measurements. 
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FIG. 1: The conditional probability distribution 
P{M; L, z) relating the galaxy luminosity L and halo mass 
M, at different redshifts, as calculated in Refs. [IB, UHl for 
SDSS at z ~ 0.05, DEEP2 at z ~ 1, and LBGs at z - 4. 
The probabilities to find a galaxy at a given luminosity in a 
halo of mass M at redshift z is plotted as a function of the 
halo mass for luminosity values for which we have galaxy 
bias data. In each of the distributions, the peak at low 
halo masses is related to galaxies of the given luminosity 
that appear as central galaxies, while the tail extending 
to higher masses is for galaxies that appear as satellites 
in more massive halos. The width of the central peak is 
related to the scatter in the relation between luminosity 
of central galaxies and halo mass and cannot simply be 
described by a delta function relating a one-to-one corre- 
spondence between mass and luminosity |19i. i20n2l|. 



These conditional probability distributions are plot- 
ted in fig [T] as a function of redshift and luminosity. 
The distribution functions account for both central 
galaxies and the satellite galaxies, following the con- 
ditional luminosity function approach of Refs. [T^.[28|. 
Due to scatter in the relation between luminosity of 



FIG. 2: The galaxy bias- luminosity data set used in our 
analysis for the three average redshifts of SDSS {z ^ 0.05), 
DEEP2 {z ~ 1), and LBG {z ~ 3.8) in comparison 
with the bias prediction calculated for the best fit ACDM 
model. The x-axis magnitude values plotted are Mr for 
SDSS, Mb for DEEP2, and Muv for LBGs at 2 ~ 3.8. 



III. ANALYSIS 

To fit bias measurements together with CMB and 
SDSS we performed a Monte Carlo Markov Chain 
analysis using a modified version of the publicly avail- 
able code cosmoMC [38^ , with a convergence diagnos- 
tic based on the Gelman and Rubin statistic [s^ (also 
known as "R-I" statistic, where R is defined as the pa- 
rameter R is defined as the ratio between the variance 
of chain means and the mean of variances). Our cos- 
moMC runs consist of 2-3 chains typically with 15000- 
20000 points and we have for our chains R—l< O.OI, 
ignoring first 50% of the chains. 

For each cosmological model we repeated the proce- 
dure described in the previous Section to calculate the 
theoretical hg{L) relation and implemented, at each 
redshift, a relation similar to that used in Ref. to 
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FIG. 3: constraints on the parameters of the ACDM+m„ model from WMAP alone (blue), WMAP+LRG+bias data 
set at z=0.05 (green) and WMAP+LRG+all bias data sets (red). 



compare with data: 

2 _ {bth,i — b^{bdatalb*)i)'^ , , 

where bth.i is the predicted bias at given luminosity 
for a given cosmological model, bdata/b^ is the ob- 
served bias at the same luminosity with error ct;,/;,, ., 
and Usys = 0.03 is a systematic uncertainty in the 
modeling of bias (40| . The sum is over the number 
of bias data points at each redshift. &*(2:) is the bias 
parameter, as a function of redshift, that accounts for 
an overall uncertainty in the bias measurements or 
modeling of bias based on CLFs. We treat it as a 
free parameter and marginalize over it when quoting 
cosmological parameter errors. This parameter shifts 
the model (or data) by a constant factor while keep- 
ing the shape the same. Thus, cosmology is measured 
through the shape of the bg{L) relation and not from 
its exact amplitude. 

We included equati on (Ht into the likelihood for the 
five-year WMAP data 'i^nd SDSS LRG power spec- 
trum (43j . We sample first the following simple seven- 
parameters cosmological model assuming flat priors 
on parameters and treating the dark energy compo- 
nent as a cosmological constant: the physical baryon 
and cold dark matter densities, Q.bh^ and flch^, the 



ratio of sound horizon to the angular diameter dis- 
tance at decoupling, 9s, the overall normalization of 
the spectrum at A; = 0.002 h Mpc~^, As, the ampli- 
tude of SZ spectrum, Asz, the optical depth to reion- 
ization r, and the scalar spectral index Us- In our 
analysis we always assume spatial flatness {flk — 0). 
The bias parameter depends on redshift, so when us- 
ing all redshift data set we introduce three free bias 
parameters for b^ in our analysis in addition to the 
cosmological parameters listed above, for a total of 
ten free parameters. 

We also explored a larger set of parameters, intro- 
ducing the sum of neutrino masses ^ rrii, and the 
dark energy equation of state w. When including 
both neutrino masses and w we performed an anal- 
ysis combining bias data with WMAP and LRG only 
and one also using Arcminute Cosmology Bolometer 
Array Receiver (ACBAR) da ta j44l| and luminosity 
distance SN-Ia data (SNe) [45l.l46l| assuming the prior 
from Hubble Space Telescope (HST) on the value of 
Hubble constant h = 0.72 ± 0.07 [13 . 

We use galaxy bias data at three different redshifts: 
from Sloan Digital Sky Survey [3l| , six points between 
redshift 0.05 and 0.1, from DEEP2 redshift survey [H 
four points at z 0.8 to z ~ 1.1, and from clustering 
of Lyman Break Galaxies in the Subaru Deep Field 
(27I with three points at z ~ 3.8. These data 
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FIG. 4: joint two-dimensional posterior probability con- 
tour plot in the crg-ria plane showing 68% and 95% con- 
tours from WMAP alone (red) and WMAP+LRG+bias 
data at all redshifts (green). 
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FIG. 6: joint two-dimensional posterior probability con- 
tour plot in the ^ m^-w plane showing 68% and 95% con- 
tours from WMAP-f LRG and bias data at all redshifts 
(red) and WMAP-|-ACBAR-HSNe-FLRG-FHST and bias 
data at all redshifts (green). 
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FIG. 5: joint two-dimensional posterior probability con- 
tour plot in the o"8-X^ rn,^ plane showing 68% and 95% 
contours from WMAP+LRG-f bias data at all redshifts. 



points are shown in fig. [2] with the best fit model. 



IV. RESULTS 

The constraints on cosmological parameters are 
shown in the Tables HI HIl and HVl with a comparison to 
constraints from WMAP five-year data alone both 
for a simple ACDM model (Table [Til and for a model 
with a non-zero mass for neutrinos (ACDM-t-miy) (Ta- 
ble [l|, and with a dark energy equation of state differ- 
ent from the cosmological constant value of —1 in ad- 
dition to neutrino mass {ACDM+m^+w) (Table ITV)) . 



In addition to WMAP data, we also consider the com- 
bination of WMAP data and SDSS LRG power spec- 
trum shape with bg{L) relation from SDSS, and finally 
the same data complemented with high-redshift bg(L) 
relations from DEEP2 and Subaru LBGs. 

As we have discussed in the introduction galaxy 
bias depends on rms fluctuation a{M) in spheres that 
contains a mass M. Galaxy bias measurements, as 
a function of luminosity, are therefore able to con- 
strain all cosmological parameters that affect this 
quantity, mainly the amplitude of matter fluctuations 
as, power spectrum spectral index or tilt n^, and neu- 
trino mass, that affect the growth of density pertur- 
bations. As presented in Ref. [2^, bias data at low 
redshifts from SDSS are already strongly sensitive to 
neutrino masses: with WMAP first-year data com- 
bined with SDSS galaxy power spectrum shape and 
SDSS bg{L) data lead to J] < 0.54 eV at the 95% 
confidence level (See Table IIIII for recent results on 
neutrino masses). 

With WMAP 5-year data and SDSS LRG power 
spectrum complemented by galaxy bias data at z 
0.05, 1, and 3.8, we are able to improve constraints 
on the sum of neutrino mass by a factor ~ 2 with re- 
spect to the result of Ref. [1^ obtaining ^ rrii, < 0.28 
eV at the 95% confidence level. We get a similar re- 
sult if we only keep to z ~ 0.05 SDSS bg{L) data and 
the LRG power spectrum shape with WMAP 5-year 
data, since by adding additional bias data at higher 
redshifts we are also introducing to the analysis two 
more unknown parameters, i.e. the nuisance bias nor- 
malization parameters for redshifts z = I and z = 3.8, 
which are marginalized over when quoting parameter 
errors. The relative increase of a factor of ~ 2 in 
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ACDM+m 

WMAP5 


±SDSS±all z's 


ACDM±TO ± m^, 

WMAP5 


±SDSS±all z's 




0.02273 ± 0.00065 


0.02248 ± 0.00060 


0.02222 ± 0.00063 


0.02234 ± 0.00060 




0.1102 ± 0.0065 


0.1160 ± 0.0041 


0.1119 ±0.0064 


0.1223 ±0.0065 


T 


0.086 ±0.017 


0.087 ±0.017 


0.083 ± 0.016 


0.083 ±0.016 


Us 


0.963 ±0.016 


0.965 ±0.014 


0.948 ±0.017 


0.954 ±0.014 


W 


-1.06 ±0.41 


-1.12 ±0.10 


-1.23 ±0.55 


-1.30 ±0.19 


Zn(10"'A«) 


3.18 ±0.05 


3.22 ± 0.04 


3.22 ±0.06 


3.24 ±0.05 




0.27 ±0.10 


0.259 ±0.020 


0.31 ±0.13 


0.267 ±0.027 


erg 


0.81 ±0.14 


0.802 ±0.037 


0.71 ±0.14 


0.775 ± 0.045 








< 1.5eV (95%CL) 


< 0.59eV (95%CL) 






1.01 ±0.03 




1.03 ±0.04 


hi 




1.18 ±0.06 




1.21 ±0.07 


hi 




3.44 ± 0.35 




3.44 ± 0.35 



TABLE I: Mean values and Icr constraints on cosmological parameters from WMAP±SDSS+bias data at all redshifts 
in comparison with constraints from WMAP alone, for models with dark energy equation of state allowed to vary. 



the neutrino mass limit compared to Ref. 23] is part 
due to the improvement in both the CMB (WMAP 
one-year to WMAP five-year) and galaxy power spec- 
trum shape data (SDSS DR2 power spectrum with 
- 200, 000 galaxies to SDSS DR4 LRG power spec- 
trum with ^ 400, 000 galaxies) ad part due to the 
improvement in the CLF modeling of the P{M\L,z) 
relation for SDSS galaxies. While the combination of 
all bias data at the three redshifts does not improve 
the limit on the sum of neutrino masses compared 
to the case with bias measurements from SDSS only, 
we do find small improvements in the uncertainties of 
the other parameters, as shown in fig. [3] for the case of 
the ACDM model with a non-zero mass for neutrinos. 
We plot also probability contours in a^-ns plane and 
in ^ TOy-(Tg plane in figs. 3] and [51 

As the growth of structure depends also on the dark 
energy density and equation of state we explored a 
more general parameter space, relaxing the assump- 
tion of a cosmological constant for dark energy and 
constraining the equation of state of dark energy w 
both in the case of neutrino mass fixed to zero and al- 
lowed to vary (Table|T|. We find that w = -1.06±0.41 
with WMAP 5-year data alone and -1.12 ± 0.10 with 
WMAP 5-year+SDSS LRG power spectrum shape 
and all hg{L) data. For comparison, the WMAP 5- 
year data combined with Baryon Acoustic Oscillation 
(BAG) data gives w = -1.15 ± 0.21. 

Our constraints on erg are erg = 0.759 ± 0.025 in the 
case KCDM+m^, while WMAP combined BAOs [6l| 
and SNe data gives erg = 0.732 ± 0.062. The analysis 
of bias data combined with CMB and SDSS showed 
in m gives 0.854 ± 0.062. 

For the most general parameter space explored in 



our analysis, with both sum of the neutrino masses 
and dark energy equation of state allowed to vary, 
we improve constraints on equation of state with 
respect to WMAP alone [w = -1.23 ± 0.54), ob- 
taining w = —1.30 ±0.19. When w is allowed to 
vary, constraints on neutrino masses are weakened to 
rrii, < 0.59 eV at the 95% confidence level, but are 
still improved with respect to ^ < 0.66eV at the 
same 95% confidence level with WMAP+BAO+SNe, 
for the same model. We then performed the same 
analysis combining WMAP, LRG and all bias data 
set with ACBAR and SNe data, assuming the prior 
from Hubble Space Telescope (HST) on the value of 
Hubble constant and we improve the last constraints 
on neutrino masses and w, obtaining ^ < 0.56eV 
at the 95% confidence level and w = -1.125 ± 0.092. 
These results are shown in table IIVI In fig [5] we plot 
the probability contours for w and ^ rrii, from WMAP 
combined with LRG and bias data sets and the same 
combined also with ACBAR and SNe, with the prior 
from HST. In the last case the inclusions of SNe and 
ACBAR data and of the prior on h slightly improves 
the constraints on ^ but above all improves con- 
straints on w (see also tables lll lIV|) . 

It is interesting to note that these limits on the sum 
of neutrino masses are almost the strongest cosmolog- 
ical constraints available in literature (as we can see 
from Table HIT)) . Compared to the results we present 
with a combination of WMAP 5-year data, SDSS 
LRG power spectrum shape, and bg{L) data at three 
redshifts, better constraints on the sum of neutrino 
masses have been published with cosmological analy- 
ses that also make use of clustering measurements of 
the Lyman-a forest [i^ [5l[ . There is some possibility 
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ACDM 

WMAP5 


+SDSS+all z's 


WMAP5 


±SDSS±all z's 




0.02273 ± 0.00062 


0.02266 ± 0.00057 


0.02226 ± 0.00063 


0.02267 ± 0.00058 




0.1099 ±0.0062 


0.1131 ±0.0034 


0.1110 ±0.0062 


0.1141 ±0.0038 


T 


0.087 ±0.017 


0.092 ±0.017 


0.084 ±0.016 


0.094 ± 0.016 


Us 


0.963 ± 0.014 


0.964 ±0.013 


0.950 ±0.017 


0.964 ± 0.013 




3.18 ±0.05 


3.20 ±0.04 


3.21 ±0.05 


3.21 ±0.04 




0.258 ±0.030 


0.273 ±0.017 


0.331 ± 0.066 


0.282 ± 0.023 




0.796 ± 0.036 


0.807 ±0.021 


0.675 ± 0.084 


0.759 ± 0.025 


Em,!/ 






< 1.3eV (95%CL) 


< 0.28eV (95%CL) 


6^ 




1.01 ±0.03 




1.03 ±0.03 


hi 




1.21 ±0.06 




1.24 ±0.06 


hi 




3.44 ± 0.34 




3.41 ±0.33 



TABLE II: mean values and Icr constraints on cosmological parameters from WMAP±SDSS±bias data at all redshifts 
in comparison with constraints from WMAP alone, for A CDM and ACDM±mi/ models. 



Kahniashvili et al. (2005) [53] Cluster Number Density 




< 


2.4 eV 


Komatsu et al. (2008) [42] 


WMAP5 




< 


1.3 eV 


Tegmark et al. (2006) [43] 


SDSS±WMAP3 




< 


0.9 eV 


Komatsu et al. (2008) [42] 


WMAP5±BAO±SNe 




< 


0.61 eV 


Kristiansen et al. (2007) [50] 


WMAP3±SDSS±SNLS±BAO±CMF 




< 


0.56 eV 


Seljak et al. (2005) [23] 


WMAP1±SDSS±SDSS hg{L) 




< 


0.54 eV 


Mac Tavish et al. (2005) [52] 


CMB±LSS 




< 


0.48 eV 


Seljak et al. (2004) [24] 


WMAPl±SDSS±Lya 




< 


0.42 eV 


Kristiansen et al. (2007) [50] 


WMAP3±SDSS±SNLS±BAO±HST 




< 


0.40 eV 


Fogli et al. (2008) [49] 


CMB±HST±SNe±BAO±Ly-Q 




< 


0.19 eV 


Seljak et al. (2006) [51] 


CMB±SDSS±2dF±SNe±Ly-a 




< 


0.17 eV 


This paper 


WMAP5±SDSS LRG± SDSS, DEEP2, LBG bg{L) Y^m^ 


< 0.28 eV 



TABLE III: Summary of the constraints at 95% confidence level on the sum of neutrino masses from various data sets 
in the literature. CMB means the collection of CMB data sets listed respectively in [49l |. [sij and [H^l. LSS means 
combination of SDSS and 2dF data sets. SNLS is the Supernova Legacy Survey, CMF is the cluster mass function, and 
Ly-Q are the clustering measurements of the Lyman-a forest. We refer the reader to individual references given in the 
table for additional details related to datasets used and how constraints were dervied. 



that Ly-a statistics may be more subject to uncertain- 
ties in both the measurement and the modehng from 
the theory side. 

The strong limit imposed with Ly-a measurements, 
combined with other cosmological data, that lead to 
^m,^ < 0.17 eV (Sll, rules out the evidence for a 



non-zero neutrino mass claimed in Ref. |58| , where the 
combination of CMB data with 2 degree Field Galaxy 
Redshift Survey [s^ data, X-ray luminosity function 
observations (XLF) and baryonic gas mass fraction 
measurements led to the constrain J2 ''^v = 0-56to;26 
eV. Independent of Ly-a data, with a 95% confidence 
level limit on the sum of neutrino masses of 0.28 eV, 
we can also conclusively state that this suggested de- 
tection of neutrino masses limit is not compatible with 



a different set of cosmological data involving galaxy 
clustering and clustering bias, combined with CMB. 
The origin of the result in Ref. [58[ was due to the 
inclusion of X-ray measurements that favored low val- 
ues of matter amplitude fluctuations, cts = 0.70±0.04. 
Due to the degeneracy between ^ and (Tg from 
CMB data the inclusion of XLF data in the analy- 
sis led therefore to the evidence for a nonzero neu- 
trino mass, that can be ruled out when considering 
increased uncertainties in the X-ray data. We also 



refer the reader to Refs. 
discussions about this result 



and i60] for further 



While our limit on 0.28 eV, in the near future, cos- 
mological data could reach the sensitivity level of ~ 
0.1 eV to begin to distinguish between the normal and 
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the inverted neutrino mass hierarchy. For example, 
inverted hierarchy could be ruled out if we can ex- 
clude that ^niy > 2y^A 77123) where A 7*7.23 ~ 2.4 • 
lQ~^eV'^ is the squared mass difference between neu- 
trino mass eigenstates [s^ . We also note that the cur- 
rent limit on the sum of neutrino masses from cosmo- 
logical observations could be quite important for ongo- 
ing and future experiments which aim to measure neu- 
trino masses: for example there is a tension between 
our results and the limits of 0.16 < miip < 0.52(2(t) 
on the neutrino mass coming from the analysis of part 
of the Heidelberg- Moscow experiment on neutrinoless 
double beta decay [56| : moreover, the next-generation 
tritium /?-decay experiment KATRIN f57] would not 
be able to measure the absolute value of neutrino mass 
because its detectability threshold is at ^ 0.2 eV. 



WMAP5-HSDSS 
+ACBAR+HST-|-SNe-hall z's 





0.02264 ± 0.00054 




0.1203 ±0.0055 


r 


0.088 ±0.016 


ris 


0.961 ±0.13 


w 


-1.125 ±0.092 




3.22 ±0.05 


O771 


0.285 ± 0.023 


o-g 


0.778 ± 0.037 




< 0.56eV (95%CL) 


bl 


1.03 ±0.02 


bl 


1.22 ±0.07 


bl 


3.41 ±0.33 



TABLE IV; Mean values and la constraints on cosmo- 
logical parameters from WMAP combined with ACBAR 
experiment, SNe — la data, HST prior and all bias data 
sets. 



V. CONCLUSIONS 

Here we have used galaxy bias measurements as a 
function of luminosity to put constraints on cosmo- 
logical parameters, mainly cts, ^tm,^, and the dark 
energy equation-of-state w. We used three galaxy 
bias-luminosity data sets at three different redshifts 
and improved by roughly a factor of two previous con- 
straints obtained with an analogous set of data, find- 
ing the sum of neutrino masses to be ^ < 0.28eV 
at the 95% confidence level for a KCDM + model, 
with (78 — 0.759±0.025. We also have shown that red- 
shift evolution of the bias information can constrain 
the equation-of-state of dark energy and we obtained 
the constraints w = -1.12 ± 0.10 for a KG DM + w 
model and 77; = -1.30 ± 0.19 for a AC DM + Tn,y + w 
model. 

In the case of dark energy equation of state allowed 
to vary the constraints on the sum of neutrino masses 
are weakened to ^rrii, < 0.59eV at the 95% confi- 
dence level, but this still improves previous constraints 
from WMAP combined with BAO and SNe data only 
(X)"^!^ < 0.66 eV). While the inclusion of Ly-alpha 
data have led to stronger constraints on the sum of 
neutrino masses, our constraint with a minimal set of 
data is still competitive and only uses clustering infor- 
mation of galaxies and the primordial fluctuations as 
probed by the CMB. In the future methods such as the 
one we use could further improve constraint on neu- 
trino masses independent of Ly-a data and achieve 
the sensitivity necessary to distinguish between the 
normal and the inverted neutrino mass hierarchies. 
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